Quantum Fluctuations to Cause the Breakdown of the Spin-1 Haldane Phase 
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Investigating quantum fluctuations in the ground states of S = 1 quantum antiferromagnetic 
spin chains described by the bilinear-biquadratic Hamiltonian Tt — y\ \Si ■ Si+i + (3(Si ■ Si+i) 2 ] , 
we study a mechanism of the breakdown of the Haldane phase. Based on the valence-bond-solid 
structure, but replacing two links of them by triplet bonds (crackions) , we construct a trial wave 
function which is singlet and translationally invariant, where the crackion-crackion distance is re- 
garded as a variational parameter. At j3 < 1/3, the minimization of the variational energy results 
in a bound state of crackions, while at /3 > 1/3, crackions come to be set free from their bound 
state with increase of f3 and the chain length. We point out that the breakdown of the Haldane 
phase with f3 approaching 1 can be attributed to the collapse of the bound state and the growth of 
a short-range repulsive interaction between crackions. 

PACS numbers: 75.10.Jm, 75.50.Ee,75.60.Ch 
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I. INTRODUCTION 

Since Haldane |Q proposed that the one-dimensional 
spin-5 Hcisenberg antiferromagnet should exhibit quali- 
tatively different properties according to whether S is in- 
teger or half odd integer, low-temperature properties of 
integer-spin chains have been of great interest. Various 
numerical tools d]|,||,[|,U brought us precise estimates of 
the excitation gap immediately above the ground state 
and visualized the Haldane massive phase. On the other 
hand, a rigorous treatment |?J performed by Affleck, 
Kennedy, Lieb, and Tasaki (AKLT) considerably con- 
tributed to understanding of the physical mechanism of 
this phenomenon. They considered a special bilinear- 
biquadratic Hamiltonian of S = 1, 



Si ■ Si+i + ~^{Si ■ Si+i) 



(1.1) 



and constructed its ground state using valence bonds, 
which they call a valence-bond-solid (VBS) state. This 
model possesses the unique disordered ground state with 
a finite gap to the excited states and therefore exhibits 
the typical nature peculiar to the Haldane phase. Keep- 
ing in mind that within a naive variational treatment 
P,p| pO|JTl| , p^ | the ground states of a certain class of 
Hamiltonians are all approximated by the VBS state, we 
are convinced that the appearance of the Haldane phase 
is generally described in terms of the VBS picture. 

The AKLT model reveals that there lies a hidden topo- 
logical order Jl3|,[l4| in the Haldane phase, which can be 
measured by the string order parameter of den Nijs and 
Rommelse [|L5|, 



O: 



string 



lim {SfJlexplmSfiS; 



(1.2) 



where ( ) denotes the expectation value in the ground 
state. Since the string order parameter takes its full value 



4/9 in the VBS state, any excitation of the VBS state 
should more or less reduce the hidden antiferromagnetic 
order. A local defect in the VBS state is obtained by re- 
placing its arbitrary link by a triplet bond jl6| . This spin- 
1 excitation, which is called a crackion, has a solitonic na- 
ture. Actually, a moving crackion and a moving domain 
wall in the hidden antiferromagnetic order both result in 
the same dispersion relation [ |17|Jlq| . Furthermore, Scharf 
and Mikeska jl!| reported that modified crackions, which 
they call dressed solitonic excitations, almost perfectly 
reproduce the low-lying excitations of the AKLT model 
Jl8| . Not only variational treatment s |2(]^j] , ^2| but also 
numerical investigations P , |2l||2^ , p5| , [2 6|]27| showed that 
the idea of domain walls appearing in the hidden order 
well applies to the low-lying excitations of the Heisenberg 
Hamiltonian as well. 



Thus the low-energy structure of S = 1 antiferromag- 
netic spin chains is well described by the VBS ground 
state and the solitonic excitations in its hidden topo- 
logical order. However, we have, on the other hand, to 
recognize that once the Hamiltonian deviates from the 
AKLT point, quantum fluctuations begin to destroy the 
perfect hidden order. At the Heisenberg point, for ex- 
ample, O s z tring is reduced to 0.374 @f|£7). Although 
the VBS model elucidates the essential mechanism of 
the appearance of the Haldane phase, the breakdown of 
the Haldane phase is never understood without consid- 
eration of the quantum fluctuations. Since a density- 
matrix renormalization-group calculation |27| suggests 
that 0| tr ; ng persists as far as the system lies in the Hal- 
dane phase, the breakdown mechanism must be revealed 
in clarifying how the string order is reduced and disap- 
pears. Why is 0| tring able to remain finite in spite of 
fluctuations? What is the driving force of the collapse of 
the string order? In answering these questions, we have 
first to recognize that there exist two mechanisms to de- 
stroy the string order. One may be called the thermal 
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mechanism and explains how the string order vanishes at 
finite temperatures. The other may generically be called 
the quantum mechanism and explains how the ground- 
state phase transitions occur. The qualitative difference 
between them was pointed out by the present author and 
Miyashita Making use of quantum Monte Carlo 

snapshots, they demonstrated that thermal fluctuations 
can cause the collapse of the string order, whereas quan- 
tum fluctuations no more than reduce the string order 
(Fig. [I]). It is worth mentioning that thermal and quan- 
tum fluctuations should be distinguished only in integer- 
spin chains. Figures 0(c) and 0(d) show that for 5=1/2 
both thermal and quantum fluctuations merely result 
in a pair of domain walls in the antiferromagnetic or- 
der and therefore there is no qualitative difference be- 
tween them. Thus we are well convinced that thermal 
fluctuations induced in the 5 = 1 Haldane phase in- 
deed lead to the collapse of the string order. On the 
other hand, quantum Monte Carlo observations [Q of 
the ground state of the Heisenberg Hamiltonian suggest 
that there exists an attractive interaction between topo- 
logical defects, which confines quantum fluctuations to 
producing local effects on the hidden order. Furthermore 
the present author [^8| analytically demonstrated that 
quantum-machanically induced solitonic excitations are 
stabilized into a strongly bound state at the Heisenberg 
point. That is why the string order parameter persists, 
in spite of fluctuations, in the ground state of the 5=1 
Heisenberg model. 

Now we have to recognize that the quantum break- 
down of the Haldane phase has remained unexplained. 
How does the bound state of the solitonic excitations be- 
have as the Hamiltonian approaches critical points? This 
is the interest in this article. Because of two reasons, we 
take the 5 = 1 bilinear-biquadratic Hamiltonian 



n = Y,[s t -s x+1 + p(s t -s l+1 ) 2 ] 



(1.3) 



on a ring of L sites for our subject. First, the model 
contains the AKLT point, which allows us to make a 
variational approach at the idea of crackions appearing 
in the VBS background. Second, the model contains ex- 
plicit critical points, which were revealed by the Bethe- 
ansatz method. One is the Uimin-Lai-Sutherland (ULS) 
point [||j3^j3ltl of = 1 and the the other is Takhtajan- 
Babujian (TB) point p2| , |3"3| ] of (3 = — 1. Hence we can 
take advantage of the knowledge on the boundaries of 
the Haldane phase in investigating the mechanism of its 
breakdown. It is the growth of quantum fluctuations with 
(3 moving from 1/3 to 1 that we explicitly discuss here. 
Another boundary of the Haldane phase, = — 1, seems 
to be so far from the AKLT point as not to allow us 
to approach it with the same scenario as demonstrated 
in this article. We will briefly mention another possible 
scenario to destroy the string order in the final section. 

Although quantum Monte Carlo snapshots are quite 
helpful in getting a qualitative view of fluctuations p4| , 
they are not naively available in the region of [3 > 
due to the negative-sign problem. Therefore, in order 



to investigate quantum fluctuations, we make an ana- 
lytic approach constructing a physically-motivated trial 
wave function for the ground state. Here, we take lit- 
tle interest in obtaining a superior variational bound jicf l 
on the ground-state energy but lay a great emphasis on 
clarifying how the pair-crackion fluctuations grow as f3 
moves away from 1/3. The trial wave function is naive 
but interestingly suggests a probable mechanism for the 
breakdown of the Haldane phase caused by quantum fluc- 
tuations. 



II. TRIAL WAVE FUNCTION 

Before constructing a trial wave function, we briefly 
review our knowledge on the ground state of the model 
as a function of (3. For recent years, the model has vig- 
orously been argued and up to now turned out to exhibit 
at least three different phases: 



— 1 < [3 < 1: Haldane phase |1|,|7J with a unique dis- 
ordered ground state and a gapped spectrum. 

(3 < —1: Dimerized phase 34 3^j3(| with twofold 
degenerate ground states and a gapped spec- 
trum. 



1 < (3: Trimerized phase |? 37 3^,||j| with threefold 
degenerate ground states. Whether the spec- 
trum is gapped or gapless has less been settled 
so far. 

A qualitative description of the phase diagram is ob- 
tained by a simple variational treatment j|,|lO|,[38| em- 
ploying three different balence-bond states, 

|VBS(L)) = Tr \g{ ® «?f ® • • • ® g s L ] , (2.1a) 
|Dimcr(L)) = Tr [f[ ® /' ® • • • ® /£] , (2.1b) 
|Trimer(L)) = Tr [h{ ® h s 2 ® • • • ® h" L ] , (2.1c) 

where 

-10), -V2|+) 
V2\-)i 10), 

= [ |->2i-i |0> 2i -i |+)ai-i ] , 

|+>* -|0)m I— >a* ] , 

V2|+) 3 ^2 |0)3i_ a v / 2|-) 3 »-2 ] , 

v / 2|-) 3 i-i -|0)3i_l 

-v / 2|-) 3i -i v / 2|+) 3l -i 

|0>3i-i -V2|+) 3 *-i 

h s 3i = T [ V2\+) 3i |0) 3l v / 2|-) 3 i ] , 



91 = 

fii-l = 

/I, - T 

h s — 
n 3i-2 — 



(2.2a) 



(2.2b) 



(2.2c) 



with |0)j, being the 5f eigenstates for eigen- 

values 1, 0, -1, respectively. |VBS(L)) is the VBS 
state of AKLT type, which is homogeneous. |Dimer(L)) 
and |Trimer(L)) represent balence-bond states which are 
dimerized and trimerized, respectively. The linear com- 
bination of these three states, 



|3(L;M)> 



cos ( 



|VBS(L)) 
|VBS(i)| 



sin (7 cos ( 



|Dimcr(i)) 
||Dimer(L)|| 



sin # sm < 



|Trimer(L)) 
||Trimer(L)|| 



(2.3) 



may be a naive but suggestive variational wave function for the ground state of the present Hamiltonian, where \\A\\ 
denotes the norm of the state vector \A). Because of the asymptotic orthogonality between the different balence-bond 
states, \<f>(L;0, (j))) is correctly normalized in the thermodynamic limit. The variational energy is obtained as 



lim 



($(L;M)|«|$(L;M)> 



4 

-3+2^- 



(2.4) 



which leads to a simple solution, 



7T 
1 



7T 
"2" 



■ = 



09 < -J). 



(|</3)- 



(2.5) 



Thus we expect that in a certain region around the 
AKLT point, the low-energy physics of the model may 
be described in terms of the VBS picture. Although 
the model encounters a commensurate-incommensurate 



crossover |27]|I(J at a certain positive value of [3, the 
gapped spectrum and the string order parameter both 
persist in the whole region between the TB (j3 = — 1) 
and the ULS {fi = 1) points fF§. 

We propose an idea of recognizing the breakdown of 
the Haldane phase with (3 moving away from the AKLT 
point as the growth of quantum fluctuations in the hid- 
den antiferromagnetic order. The crackion with its spin 
projection A (A = +, 0, — ) is created at the bond between 
sites i and i + 1 replacing gf by g£ Jill , where 



9 = 



V2\+) 
-|0) 



9° = 



-|0) V2\+) 
V2\-) -|0) 



-|0) 
V2\-) 



(2.6) 



The site indices have been omitted in Eq. (2J3) for the sake of simplicity. In order to describe the ground state, we 
introduce a trial wave function of zero momentum with a pair of crackions, 



\^(L;l,a)) = ^=Yl 



I* 



(£)> 



I* 



(L)) 



I* 00 , 



11^^)1 



(L)\\ 



(L)) 



\\*%+i(L)\\ 



(2.7) 



where 



l*fc(£)> = 



Tr [g\ ® • • • ® g!_ ± ® 5j A ® gf +1 ® • • • 

for 1 < I < L/2, 
Tr [gf <g> • • • ® 5 fJ for Z = . 



(2.8) 

Here, the crackion-crackion distance Z is regarded as a 
variational parameter, whereas a has just been intro- 
duced so as to symmetrize the two crackions into a spin 
singlet. The condition of the total spin singlet, 



L 

results in 



V(L;l,a)\Si ■ Sj\^(L;l,a)) = 0, (2.9) 



a = aa 



y 3 L + (-l) l 3 L - 1 + {-l) L - l 3 l + (-1) L ' 

(2.10) 

We note that «o generally deviates from unity. This is 
because a crackion has its internal structure and spreads 
over neighboring sites. It is convincing that «o co- 
incides with unity only in the limit of L — ► oo and 



Z(< L/2) — > oo. On the other hand, the minimization 
of the variational energy with respect to a, 

d (y(L;l,a)\H\V{L;l,a)) 



da 



0. 



(2.11) 



(<I>(L;l,a)\y(L;l,a)) 
also results in Eq. (2.10). In this sense the trial wave 



function employed is reasonable enough to investigate the 
ground-state fluctuations of the present Hamiltonian. 

Hereafter we fix a to ctQ and calculate the variational 
energy as a function of L, Z, and /3, 

{*(L-,l,ao)\HMm,ao)) 



E(L;l,/3) 



(2.12) 



(^(L;l,a )\^(L;l,a )) 

The stabilization of the two-crackion state is measured 
by the energy difference 

AE(L; Z, 0) ee E(L; Z, (3) - E(L; 0, (3) . (2.13) 



quantum averages ('J^'f. i (L) \ S 



The variational energy (2.12) is constructed from the 



A f t. 



S m+1 \*% +l (L)) 



and 



(L)\(Sr. 



U,i+l\^JW~"m ■ S m+ i) 2 \^™ +l (L)), which are repre- 
sented in the form of matrix-product type and therefore 
calculated by the use of the transfer-matrix technique 
p8| , |41 42 ,4^]. The calculation is straightforward but ex- 
hausting due to plenty of indices appearing in the matrix 
elements. The explicit form of the variational energy is 
too lengthy to be presented even though a is fixed. 
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III. RESULTS 



tion of the ground state: 



First, let us concentrate our attention on the Heisen- 
berg point (3 = 0. We list in Table | the stabiliza- 
tion energy of the two-crackion state, AE(L; 1,0). With 
AE{L\ I, 0) < 0, a pair of crackions with their distance 
I may spontaneously occur in a chain of length L. The 
crackion fluctuation is even unstable in too short chains. 
Actually, at L < 8, the VBS state itself turns out to be 
the variational bound. However, a pair of crackions come 
to be stabilized into a strongly bound state with increase 
of L. The extremely small energy stabilization at I = 1 
may be attributed to the fact that a crackion spreads 
over two sites. As L increases, the crackion pairs with 
larger values of I are stabilized one after another. How- 
ever, the maximum distance between crackions keeping 
their formation energy negative stays four at L = 1000 
and five even at L = 10000. Therefore two crackions can 
hardly move away from each other. Regarding L, as well 
as I , as a variational parameter, we obtain the variational 
bound on the ground-state energy per site, —1.37012987, 
at L — 56 and 1 — 2. Further increase of the chain length 
only reduces the crackion effect and ends in no stabiliza- 



lim 



E(L;l) 

z 



lim 



(VBS(L)|H|VBS(L)) 
L(VBS(L)|VBS(L)) 



The ground state in the long-chain limit should be de- 
scribed by a variational wave function with a finite 
crackion density. Alternatively, within an approxima- 
tion neglecting the interaction between the bound crack- 
ion pairs, we may conclude that the best density of the 
crackion pairs is about 1/56. Thus we are led to a likely 
ground-state picture that pairs of crackions with I ~ 2 ap- 
pear in the VBS background keeping their density about 
1/56. The true ground state is supposed to fluctuate 
around this picture. In fact the present variational bound 
has not yet reached the correct value —1.402(1), which is 
a quantum Monte Carlo estimate of the per-site ground- 
state energy at I = 56. The rest of the correlation 
may be attributed to a spatial extension of a crackion, 
a breathing motion of two crackions within a pair, and 
an interaction between crackion pairs. A crackion is not 
equivalent to a hidden domain wall but indeed possesses 
a solitonic nature For example, there is a relation 

between a pair of hidden domain walls and crackion pairs: 



Sf Sf +l \VBS(L)) = -|*£.WiW> 



(L)} 



(3.2) 



where we assume that I > 2. Considering that spin flips 
in pairs keeping the total magnetization constant do not 
necessarily cause a soliton-antisoliton pair in the hidden 
order (Fig. |2|), the present variational calculation is well 
consistent with a direct observation of the ground state 
(Fig- 1). 

Next, let us turn on the biquadratic interaction be- 
tween neighboring spins. We plot in Fig. ^ the sta- 
bilization energy of the two-crackion state, AE(L;l, (3), 
as a function of I changing L and (3. The singularity 
observed around L = 1/2 should be attributed to the 
periodic boundary condition. We find that the bound 
state of crackions, which is fully stabilized at the Heisen- 
berg point, becomes less stable with increase of (3. It is 
a matter of course that any crackion fluctuation causes 
an increase of the energy at the AKLT point (3 = 1/3, 
where the VBS state with the perfect hidden order is the 
exact ground state. Interestingly, in long enough chains 
at (3 > 1/3, a pair of crackions are less stable in a bound 
state than when they are far away from each other. As (3 
moves away from 1 /3 toward 1 as well as with increase of 
L, the attractive interaction between crackions changes 
into a repulsive one at small values of /. Although the 
stabilization energy per site is more and more reduced as 
L increases due to the reduction of the crackion density, 
yet the biquadratic interaction indeed turns the forma- 
tion energy of well-separated two crackions negative in 
long enough chains. Hence a pair of crackions are set 
free from their bound state and are stabilized in their 
wide-range breathing motion with large enough values of 



f3 in long enough chains, as is illustrated in Fig. pi We 
have already stated that the strong attractive interaction 
between crackions could allow the string order parameter 
to remain finite against quantum fluctuations. The col- 
lapse of the bound state of crackions and the growth of 
the short-range repulsive interaction between crackions, 
this is a probable mechanism of the quantum mechanical 
breakdown of the Haldane phase. 



IV. SUMMARY AND DISCUSSION 

We have investigated quantum fluctuations in the 
ground state of the spin-1 chain with the bilinear- 
biquadratic Hamiltonian. With the help of quantum 
Monte Carlo observations of the ground state, we have 
constructed a trial wave function with a pair of crack- 
ions appearing in the VBS background. We have demon- 
strated that two crackions are indeed stabilized into a 
strongly bound state at the Heisenberg point. While 
Monte Carlo snapshots are not available in the region of 
[3 > due to the negative-sign problem, we have analyt- 
ically revealed the scenario for the quantum mechanical 
collapse of the string order with [3 approaching 1. The 
disappearance of the attractive interaction sets crackions 
free from their bound state and the growth of a repul- 
sive interaction causes crackions to move away from each 
other. 

Besides the biquadratic interaction, there are many 
other factors to cause the breakdown of the Haldane 
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phase, such as an alternating or anisotropic interaction 
@|2 . Therefore there may be any other 

scenario for quantum phase transitions. Actually we can 
not describe the Haldane-to-dimer phase transition with 
[3 approaching —1 by the same scenario. As (3 moves 
from 1/3 to —1, the bound state of crackions remains 
fully stabilized. The TB point /3 = — 1 might be too 
far from the AKLT point (3 = 1/3 to be discussed shar- 
ing the same physics with the VBS Hamiltonian. How- 
ever, the picture of crackion pairs moving in the VBS 
background is still totally valid at the Heisenberg point 
p0|,p4|| . Even though two crackions are stabilized into a 
bound state, bound crackion pairs with a high enough 
density can cause the collapse of the string order. The 
breakdown mechanism of the Haldane phase with (3 ap- 
proaching — 1 may be beyond the present approximation 
neglecting any interaction between crackion pairs. Since 
it is hardly feasible to perform a direct calculation of 
numbers of crackions keeping a definite picture of fluc- 
tuations, we expect alternative approaches (n],^0| in the 
region of (3 < 0. Compared with the thermal break- 
down, quantum phase transitions are various and com- 
plicated. We hope the present argument will motivate 
further study on the quantum mechanical breakdown of 
the Haldane phase. 
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romagnetic order, (a) A thermal fluctuation in the S = 1 
system, which allows the total magnetization to change and 
therefore leads to the collapse of the long-range hidden anti- 
ferromagnetic order, (b) A quantum fluctuation in the S = 1 
system, which keeps the total magnetization constant and 
therefore only produces a local effect on the hidden antifer- 
romagnetic order, (c) A thermal fluctuation in the S = 1/2 
system, which is, in contrast with the S = 1 case, nothing 
more than a local effect on the antiferromagnetic order, (d) 
A quantum fluctuation in the S = 1/2 system, which locally 
disturbs the antiferromagnetic order. 

FIG. 2. Quantum fluctuations in the S = 1 ground state, 
where the total magnetization is kept constant and therefore 
spins flip in pairs. The black triangle denotes a domain wall 
in the hidden antiferromagnetic order. In the case (a), a soli- 
ton-antisoliton pair appears in the hidden order, whereas in 
the case (b), there appears no defect in the hidden order. 

FIG. 3. A quantum Monte Carlo snapshot of the trans- 
formed two-dimensional Ising system of L = 100 at the S = 1 
Heisenberg point under the temperature feT = 0.02, where 
the horizontal and the vertical axes denote the chain and the 
Trotter directions corresponding to space and time, respec- 
tively. The temperature taken is low enough to represent the 
ground-state properties. We have set the Trotter number n 
equal to 48 and show the passage of time corresponding to 
1/4£;bT. We encircle the antiphase domains which have the 
hidden antiferromagnetic order opposite to the background. 
The location of the domain walls is specified inevitably with 
some uncertainty owing to the liquid-like nature of the VBS 
state. 

FIG. 4. The stabilization energy of the two-crackion state 
as a function of the crackion-crackion distance / at various val- 
ues of the chain length L and the biquadratic interaction (5, 
AE(L;l,(3) = E(L;l,/3)-E(L;0,/3). (a) L = 80, (b) L = 160, 
(c) L = 240, (d) L = 320. 

FIG. 5. Schematic representation of the ground-state spin 
configurations at distinct regions of the biquadratic interac- 
tion (5. The black triangle denotes a domain wall in the hid- 
den antiferromagnetic order, (a) The AKLT point /3 = 1/3. 
We observe a perfect hidden antiferromagnetic order, (b) 
P < 1/3. The domain walls are stabilized into a bound pair. 
The hidden antiferromagnetic order is reduced but still per- 
sists, (c) f3 > 1/3. The domain walls are set free from their 
bound state. The hidden antiferromagnetic order is nonlo- 
cally broken. 
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TABLE I. The stabilization energy of the two-crackion state as a function of the crackion-crackion distance I and the chain 

length L at the Heisenberg point, AE(L; I, 0) = E(L; I, 0) - E(L; 0, 0). 

I L — 20 L = 40 L = 60 L = 80 L = 100 L = 120 L = 140 L = 160 L = 180 L = 200 

1 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 

2 0.20512 -1.33333 -2.20290 -2.76190 -3.15152 -3.43860 -3.65891 -3.83333 -3.97484 -4.09195 

3 2.21638 1.43089 0.73563 0.11594 -0.43986 -0.94118 -1.39564 -1.80952 -2.18803 -2.53552 

4 2.57370 2.44126 2.31010 2.18079 2.05330 1.92758 1.80360 1.68132 1.56071 1.44173 

5 2.72694 2.71098 2.69141 2.67187 2.65236 2.63289 2.61345 2.59403 2.57465 2.55530 

6 2.67604 2.69500 2.69236 2.68971 2.68706 2.68441 2.68177 2.67912 2.67648 2.67383 

7 2.60728 2.71581 2.71546 2.71511 2.71476 2.71441 2.71407 2.71372 2.71337 2.71302 

8 2.28123 2.70902 2.70898 2.70893 2.70889 2.70884 2.70880 2.70875 2.70871 2.70866 
1 L = 220 L = 240 L = 260 L = 280 L = 300 L = 320 L = 340 L = 360 L = 380 L = 400 

1 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 

2 -4.19048 -4.27451 -4.34703 -4.41026 -4.46586 -4.51515 -4.55914 -4.59864 -4.63430 -4.66667 

3 -2.85564 -3.15152 -3.42579 -3.68075 -3.91837 -4.14035 -4.34820 -4.54321 -4.72655 -4.89922 

4 1.32436 1.20856 1.09429 0.98153 0.87025 0.76042 0.65201 0.54499 0.43934 0.33503 

5 2.53599 2.51670 2.49744 2.47822 2.45903 2.43986 2.42073 2.40163 2.38256 2.36352 

6 2.67119 2.66854 2.66590 2.66325 2.66061 2.65797 2.65533 2.65269 2.65005 2.64740 

7 2.71267 2.71232 2.71197 2.71163 2.71128 2.71093 2.71058 2.71023 2.70988 2.70954 

8 2.70862 2.70858 2.70853 2.70849 2.70844 2.70840 2.70835 2.70831 2.70826 2.70822 
1 L = 420 L = 440 L = 460 L = 480 L = 500 L = 520 L = 540 L = 560 L = 580 L = 600 

1 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 

2 -4.69617 -4.72316 -4.74797 -4.77083 -4.79198 -4.81159 -4.82984 -4.84685 -4.86275 -4.87764 

3 -5.06215 -5.21612 -5.36185 -5.50000 -5.63113 -5.75578 -5.87440 -5.98742 -6.09524 -6.19820 

4 0.23203 0.13033 0.02989 -0.06931 -0.16728 -0.26406 -0.35966 -0.45411 -0.54742 -0.63962 

5 2.34451 2.32553 2.30658 2.28767 2.26878 2.24992 2.23109 2.21230 2.19353 2.17479 

6 2.64477 2.64213 2.63949 2.63685 2.63421 2.63157 2.62894 2.62630 2.62366 2.62103 

7 2.70919 2.70884 2.70849 2.70814 2.70779 2.70745 2.70710 2.70675 2.70640 2.70605 

8 2.70817 2.70813 2.70808 2.70804 2.70799 2.70795 2.70791 2.70786 2.70782 2.70777 
I L — 620 L — 640 L — 660 L — 680 L = 700 L = 720 L — 740 L — 760 L — 780 L — 800 

1 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 

2 -4.89162 -4.90476 -4.91715 -4.92884 -4.93989 -4.95035 -4.96028 -4.96970 -4.97865 -4.98718 

3 -6.29662 -6.39080 -6.48101 -6.56749 -6.65047 -6.73016 -6.80674 -6.88041 -6.95131 -7.01961 

4 -0.73072 -0.82075 -0.90973 -0.99767 -1.08459 -1.17051 -1.25544 -1.33941 -1.42242 -1.50450 

5 2.15608 2.13740 2.11876 2.10014 2.08155 2.06299 2.04446 2.02596 2.00749 1.98904 

6 2.61839 2.61576 2.61312 2.61049 2.60786 2.60522 2.60259 2.59996 2.59733 2.59470 

7 2.70570 2.70536 2.70501 2.70466 2.70431 2.70396 2.70361 2.70327 2.70292 2.70257 

8 2.70773 2.70768 2.70764 2.70759 2.70755 2.70750 2.70746 2.70741 2.70737 2.70733 
I L — 820 L — 840 L — 860 L — 880 L = 900 L = 920 L — 940 L — 960 L — 980 L = 1000 

1 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 

2 -4.99531 -5.00306 -5.01046 -5.01754 -5.02432 -5.03081 -5.03704 -5.04301 -5.04875 -5.05426 

3 -7.08544 -7.14894 -7.21022 -7.26941 -7.32660 -7.38190 -7.43540 -7.48718 -7.53733 -7.58592 

4 -1.58567 -1.66592 -1.74529 -1.82379 -1.90142 -1.97821 -2.05417 -2.12931 -2.20364 -2.27719 

5 1.97063 1.95225 1.93389 1.91556 1.89727 1.87900 1.86076 1.84255 1.82436 1.80621 

6 2.59207 2.58944 2.58681 2.58418 2.58155 2.57892 2.57629 2.57367 2.57104 2.56841 

7 2.70222 2.70187 2.70152 2.70118 2.70083 2.70048 2.70013 2.69978 2.69944 2.69909 

8 2.70728 2.70724 2.70719 2.70715 2.70710 2.70706 2.70701 2.70697 2.70692 2.70688 
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